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The elliptic flow excitation function calculated in a full (3+l)d hybrid Boltzmann approach with 
an intermediate hydrodynamic stage for heavy ion reactions from GSI-SIS to the highest CERN-SPS 
energies is discussed in the context of the experimental data. In this study, we employ a hadron gas 
equation of state to investigate the differences in the dynamics and viscosity effects. The specific 
event-by-event setup with initial conditions and freeze-out from a non-equilibrium transport model 
allows for a direct comparison between ideal fluid dynamics and transport simulations. At higher 
SPS energies, where the pure transport calculation cannot account for the high elliptic flow values, 
the smaller mean free path in the hydrodynamic evolution leads to higher elliptic flow values. In con- 
trast to previous studies within pure hydrodynamics, the more realistic initial conditions employed 
here and the inclusion of a sequential final state hadronic decoupling provides results that are in line 
with the experimental data almost over the whole energy range from _Ei a b = 2 — 160^1 GeV. Thus, 
this new approach leads to a substantially different shape of the v^je scaling curve as a function of 
(1/ SdN c h/dy) in line with the experimental data compared to previous ideal hydrodynamic calcu- 
lations. This hints to a strong influence of the initial conditions for the hydrodynamic evolution on 
the finally observed vi values, thus questioning the standard interpretation that the hydrodynamic 
limit is only reached at RHIC energies. 
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Transverse collective flow is one of the earliest pre- 
dicted observables to probe heated and compressed nu- 
clear matter P, Q ■ Elliptic flow, the anisotropy param- 
eter that quantifies the momentum space anisotropy in 
the transverse plane of the outgoing particles of a heavy 
ion reaction, is a result of the pressure gradients that 
are present in the course of the evolution. Since it is 
a self-quenching effect, it is very sensitive to the early 
stage of the collision, i. e. the initial conditions and the 
mean free path during the high energy density stage of 
the evolution. 

Hydrodynamics has been proposed many years ago as 
a tool to describe collective effects in the hot and dense 
stage of heavy ion reactions where the matter might be- 
have like a locally thermalized ideal fluid [E 0, IE IE IE HI • 
The hydrodynamic description has gained importance 
over the last few years because the high elliptic flow val- 
ues that have been observed at the Relativistic Heavy 
Ion Collider (RHIC) seem c omp atible with some ideal hy- 
drodynamic predictions [E lid ]. At bombarding energies 
below the highest RHIC energies however, the same hy- 
drodynamic calculations do overpredict the elliptic flow 
values by a large amount. This success of the hydrody- 
namic model has lead to the speculation that at RHIC 
energies and beyond the system reaches thermal equilib- 
rium so quickly that a dense partonic (and nearly perfect) 
liquid could form. In this paper we investigate how ro- 
bust this interpretation is, if proper initial conditions and 
freeze-out procedure are taken into account. This is nec- 
essary because hydrodynamic results depend strongly on 
the initial and final boundary conditions that are applied 



in the calculation. 

Previous calculations of elliptic flow in hadronic trans- 
port approaches have led to the conclusion that the 
pressure in the early stage of the collision is too low 
to reproduce the high elliptic flow values measured at 
RHIC [III EE [HI]. To get a more consistent pic- 
ture of the whole dynamics of heavy ion reactions var- 
ious so called microscopic plus macroscopic hybrid ap- 
proaches have been launched du ring the last decade 
E El EE EE EE EE El, HE El, EE Hf) . Here we use the 
same technique and employ a transport approach with an 
embedded three-dimensional ideal relativistic one fluid 
evolution for the hot and dense stage of the reaction 
based on the Ultra-relativistic Quantum Molecular Dy- 
namics (UrQMD) model 0, HI[. This approach allows 
to reduce the parameters for the initial conditions and 
provides a consistent freeze-out description and allows to 
compare the different underlying dynamics - ideal fluid 
dynamics vs. non-equilibrium transport - directly. 

This integrated Boltzmann+hydrodynamics transport 
approach is applied in this paper to simulate the dynam- 
ics of heavy ion collisions and to extract elliptic flow val- 
ues from i?i a b = 2 — 160A GeV. To mimic experimen- 
tal conditions as realistically as possible the initial con- 
ditions and the final hadronic freeze-out are calculated 
using the UrQMD approach. With this ansatz, the non- 
equilibrium dynamics in the very early stage of the col- 
lision as well as the final state interactions are properly 
taken into account on an event-by-event-basis. 

Let us shortly describe the features of the present ap- 
proach. UrQMD [2E I!?], Hl| is a hadronic transport ap- 
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proach which simulates multiple interactions of ingoing 
and newly produced particles, the excitation and frag- 
mentation of color strings and the formation and de- 
cay of hadronic resonances. The coupling between the 
UrQMD initial state and the hydrodynamical evolution 
proceeds when the two Lorentz-contracted nuclei have 
passed through each other [24[ ■ Here, the spectators con- 
tinue to propagate in the cascade and all other hadrons 
are mapped to the hydrodynamic grid. This treatment 
is especially important for non-central collisions which 
are of interest here. Event-by-event fluctuations are di- 
rectly taken into account via event wise initial conditions 
generated by the primary collisions and string fragmen- 
tations in the microscopic UrQMD model. This leads to 
non-trivial velocity and energy density distributions for 
the hydrodynamical initial conditions as will be discussed 
below. 

Starting from these initial conditions a full (3+1) di- 
mensional ideal hydrodynamic evolution is performed us- 
ing the SHASTA algorithm 0, The hydrodynamic 
evolution is stopped, if the energy density e drops be- 
low five times the ground state energy density Sq (i.e. 
~ 730MeV/fm 3 ) in all cells. This criterion corresponds 
to a T-/i£-configuration where the phase transition is 
expected (approximately T = 170 MeV at [is = 0). 
The hydrodynamic fields are then transformed to par- 
ticle degrees of freedom via the Cooper-Frye equation on 
an isochronous time t- hypersurface in the computational 
frame. The created particles proceed in their evolution in 
the hadronic cascade (UrQMD) where rescatterings and 
final decays are calculated until all interactions cease and 
the system decouples. Further we refer to this kind of 
freeze-out procedure as the isochronous freeze-out (IF). 

Alternatively an approximate iso-eigentime freeze-out 
is chosen (see 30] for details). Here, we freeze out full 
transverse slices, of thickness Az — 0.2fm, whenever all 
cells in each individual slice fulfill the freeze-out crite- 
rion. For each slice we apply the isochronous procedure 
described above. By doing this one obtains a rapidity in- 
dependent freeze-out temperature without artificial time 
dilatation effects. In the following we will refer to this 
procedure as "gradual freeze-out" (GF). A more detailed 
description of the hybrid model including parameter tests 
and results for multiplicities and spectra can be found in 

Serving as an input for the hydrodynamical calcula- 
tion the equation of state (EoS) strongly influences the 
dynamics of an expanding system. In this work we em- 
ploy a hadron gas (HG) equation of state, describing a 
non-interacting gas of free hadrons [3l| . Included here 
are all reliably known hadrons with masses up to « 2 
GeV, which is equivalent to the active degrees of freedom 
of the UrQMD model (note that this EoS does not con- 
tain any form of phase transition) . The purely hadronic 
calculation serves as a baseline calculation to explore the 
effects of the change in the underlying dynamics - pure 
transport vs. hydrodynamic calculation. 

We begin our investigation with the initial conditions 
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FIG. 1: (Color online) Energy density distribution in the x—y- 
plane for one midcentral (6 = 7 fm) Pb+Pb collision at = 
4(L4 GeV. 

for the hydrodynamical evolution. Fig. Q] shows the ini- 
tial local rest frame energy density distribution in the 
transverse plane for one single Pb+Pb collision. The 
spatial anisotropy that causes the development of elliptic 
flow is nicely observed. The distribution is not smooth 
and not symmetric in any direction. In Fig. [5] the ini- 
tial velocity distribution is shown in the transverse plane. 
The value that is plotted here is the absolute value of 
the three- velocity of the hydrodynamic cells (v — \v\) 
times the local rest frame energy density in the respec- 
tive cell. In this way, one gets rid of the numerical noise 
in the almost empty cells. The velocity distribution is 
also not symmetric and fluctuates from event to event. 
As expected the velocities are higher at the edges of the 
almond shaped overlap region in x direction. In the mid- 
dle of the system the matter is almost at rest. To see how 
these space-momentum correlations transform to observ- 
ables we introduce the elliptic flow V2- 

The second coefficient of the Fourier expansion of the 
azimuthal distribution of the emitted particles {v^) is 
called elliptic flow [12, [H, HH . V2 is defined by 



U2 EE(COS[2(0-$ RP )]) = (^-^|) , (1) 

\ Px ' Py I 

where <f> is the azimuthal angle of the particle, <&rp 
is the azimuthal angle of the reaction plane and p x and 
p y are the momenta of the particle in x- and y-direction 
respectively. 

In Fig. [3] we show the excitation function of charged 
particle elliptic flow compared to data over a wide en- 
ergy range (Fig. [3]), i.e from E\ &h = 1A GeV to ^/snn = 
200 GeV. In this figure, the data and the calculation are 
not divided by further model dependent quantities and 
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FIG. 2: (Color online) Velocity distribution in the transverse 
x — y-plane for one midcentral (6 = 7 fm) Pb+Pb collision at 
£aab = 40A GeV. The absolute value of the velocity has been 
multiplied by the energy density in the correspondig cell. 



therefore it allows for a direct comparison. The symbols 
indicate the data for charged particles from different ex- 
periments. In the SPS regime the pure transport model 
calculations are quite in line with the data, especially 
with the NA49 results. Above E lah = 160 A GeV the cal- 
culation underestimates the elliptic flow. This has been 
taken as a sign that partonic degrees of freedom become 
more important at these energies. 

The smaller mean free path in the hybrid model calcu- 
lation leads to higher elliptic flow values at higher SPS 
energies even without explicit phase transition. At lower 
energies the result is in line with the transport calcu- 
lation since the hydrodynamic evolution is very short. 
The average duration of the hydrodynamic evolution in- 
creases from ~ 3 fm at low energies to around 8 fm at 
-Blab = 40A GeV and even — 12 fm at the highest SPS en- 
ergy. Please note that the crucial observation is not only 
that there is higher elliptic flow than in the transport cal- 
culation, but that the hybrid approach shows that ideal 
hydrodynamics is less than 20 % away from the exper- 
imental data. This confirms that the initial conditions 
and the freeze-out treatment have important influence 
on the results of a hydrodynamic calculation. 

The elliptic flow of pions as a function of transverse mo- 
mentum is shown in Fig. [5] As it has been stated above 
the hydrodynamic evolution leads to higher elliptic flow 
values especially at higher p t where the pure transport 
calculation underpredicts the data. For these differen- 
tial results at midrapidity the difference between the two 
freeze-out prescriptions is less pronounced than for the in- 
tegrated results. The hybrid model calculation leads also 
for the differential elliptic flow to a reasonable agreement 
with the experimental data. 
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FIG. 3: (Color online) The energy excitation function of el- 
liptic flow of charged particles in Au+Au/Pb+Pb collisions 
in mid-central collisions (b=5-9 fm) calculated at midrapidity 
(\y\ < 0.5) within the hybrid model with isochronous freeze- 
out (black full line) and gradual freeze-out (black dashed line) 
is contrasted to the pure UrQMD transport calculation (black 
dotted line). These curves are compared to data (colored sym- 
bols) from different experiments (E895J5877, NA49, CERES, 
PHENIX, PHOBOS and STAR) [UlSl- 



Finally, we replot the V2(y/s NN ) values as a func- 
tion of particle density and scaled by the eccentricity 
of the initial state. Fig. [5] shows i^/e as a function 
of (1/ S)dN c h/dy which is assumed to be a decreasing 
quantity in the investigated energy regime in ideal hy- 
drodynamics calculations. (1/ S)dN c h/dy is the charged 
particle density at midrapidity divided by the initial state 
overlap area. This way of plotting the elliptic flow exci- 
tation function allows to compare results from different 
energies and centralities at the same time. The charged 
particle multiplicity in the overlap area is the same in a 
central low energy collision as in a peripheral high energy 
collision. The calculations within the UrQMD + hydro- 
dynamics approach have been performed for midcentral 
Pb+Pb/ Au+Au collisions with 6 = 5—9 fm. The charged 
particle multiplicities and v-i values have been calculated 
at midrapidity < 0.5). For the evaluation of the rel- 
evant initial eccentricity and the overlap area we have 
stopped the calculation of 10.000 UrQMD events at the 
time of the overlap of the nuclei. The quantities of in- 
terest have been evaluated for the participating particles 
which arc defined as the nucleons which have undergone 
at least one interaction plus the newly produced parti- 
cles (the result is insensitive to the fact if only nucleons 
are considered or not) according to the following formula 
MM: 
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FIG. 4: (Color online) Elliptic flow of pions in mid-central 
(b=5-9) Pb+Pb collisions at E lab = 40A GeV and E lab = 
16CL4 GeV. The full black line depicts the hybrid model cal- 
culation with isochronous freeze-out, the black dashed line 
the hybrid calculation employing the gradual freeze-out while 
the pure transport calculation is shown as the black dotted 
line. The colored symbols display experimental data obtained 
with different measurement methods by NA49 [3^ |. 
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FIG. 5: (Color online) 172/e as a function of (1/ S)dN c h/dy 
for different energies and centralities for Pb+Pb/ Au+Au col- 
lisions. The results from mid-central collisions (b=5-9 mi) 
calculated within the hybrid model with isochronous/gradual 
freeze-out (full line with circles and dashed line with trian- 
gles respectively) are shown. Furthermore, the hybrid model 
calculation with GF is divided by a different eccentricity 
(e2)(dashed line with squares). These curves are compared to 
data depicted by colored symbols from different experiments 
(E877, NA49 and STAR ^) for mid-central collisions. The 
green full lines correspond to the previously calculated hydro- 
dynamic limits [13]. 
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andS = 7rV<.x 2 )<y 2 ) , (2) 



where the averages are taken over particles and events 
at the same time. An alternative dcfintion is sometimes 
used 
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(3) 



where the averages are taken first over particles in one 
event and then the value for e is averaged over events. In 
this way, the events with higher particle production have 
the same weight as those with lower multiplicities. 

Since it is not obvious which way of calculating the 
initial eccentricity captures the physical picture best (e.g. 
the Glauber values which are taken from experiment) 
we show both possibilities because they lead to different 
results. 

It is remarkable that the shape of the curve is substan- 
tially changed in the hybrid model calculation compared 



to the previous calculated hydro limits [37| (shown as 
horizontal lines). A similar shape has also been obtained 
in a two-dimensional hybrid calculation with a hadronic 
afterburner but simplified initial conditions and a EoS 
including a phase transition to the QGP (see 16]). The 
present calculation with ideal one fluid hydrodynamics 
with hadronic degrees of freedom as described above is 
however able to reproduce the shape of the experimental 
data points and even the magnitude at lower energies. 
Note that at very low energies the hydrodynamic stage 
is rather short and does not influence the evolution con- 
siderably. In this regime all different setups show the 
same results and are compatible with the data. Towards 
higher energies (where pure hadronic transport calcula- 
tions have too much viscosity and underpredict the data), 
the hybrid calculation leads to higher pressure gradients 
in the early stage and therefore to suitable elliptic flow 
values. Most important here is the only moderate in- 
crease in i>2 because of the more realistic treatment via 
the inclusion of initial non-equilibrium effects and a com- 
plex shape (both in coordinate and momentum space) of 
the initial energy and baryon density distribution in ad- 
dition to our sophisticated final state freeze-out. In fact, 
for the (more physical) gradual freeze-out, the V2 values 
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are reduced by more than a factor of 2 compared to ideal 
hydrodynamics with simplified initial conditions. The 
influence of the freeze-out prescription can be observed 
by comparing the isochronous to the gradual freeze-out 
scenario. The alternative calculation of the eccentricity 
(£2) leads to higher values because €2 is in general smaller 
than e. Additional changes in magnitude may be caused 
by viscosity effects during the hot and dense stage and a 
possible phase transition which might weaken the pres- 
sure gradients in the early stage. 

We have shown that the elliptic flow values at SPS en- 
ergies are in line with an ideal hydrodynamic evolution if 
a proper initial state is used and the final freeze-out pro- 
ceeds gradually. An integrated hybrid model calculation 
with initial conditions and freeze-out from microscopic 
transport with an embedded (3+1 )d hydrodynamic evo- 
lution is able to reproduce the experimentally measured 
V2 values. This points to the fact that the treatment of 
initial conditions and freeze-out is crucial for any hydro- 



dynamic calculation. 
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